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' The nonequilibrium tunneling center model of a localized electronic level coupled to a fluctuating 

, two-state system and to two electronic reservoirs, is solved via an Anderson- Yuval-Hamann mapping 

onto a plasma of alternating positive and negative charges time-ordered along the two "Keldysh" 
contours needed to describe nonequilibrium physics. The interaction between charges depends both 
^ , on whether their time separation is small or large compared to a dephasing scale defined in terms 

■ of the chemical potential difference between the electronic reservoirs and on whether their time 

separation is larger or smaller than a decoherence scale defined in terms of the current flowing from 
one reservoir to another. A renormalization group transformation appropriate to the nonequilibrium 
problem is defined. An important feature is the presence in the model of a new coupling, essentially 
the decoherence rate, which acquires an additive renormalization similar to that of the energy 
in equilibrium problems. The method is used to study interplay between the dephasing-induced 
, formation of independent resonances tied to the two chemical potentials and the decoherence which 

i-^ ■ cuts off the scaling and leads to effectively classical long-time behavior. We determine the effect of 

^ ' departures from equilibrium on the localization-delocalization phase transition. 



X 



I. INTRODUCTION 



Understanding the nonequilibrium behavior of interacting quantum mechanical systems is one of the important open 
issues in condensed matter physics, with applications in nanoscience-'^ , the study of cold atoms in optical lattices'^, 
' ^ I nonlinear spectroscopies- and, transport at quantum critical points^. One may distinguish three classes of nonequi- 
g . librium situations: response of a system initially in an equilibrium state to a strong transient pulse, time evolution of 
^ a system from a particular initial condition, and the steady state behavior of a driven system. In this paper we shall 
I— I, be concerned with one of the simplest examples of the third class of problems: a quantum mechanical system with 
only a few degrees of freedom, coupled to two reservoirs with which particles and energy may be exchanged, and with 
^ ] the nonequilibrium drive arising from a difference, A/x, in chemical potential n between the reservoirs. This model 
I ' is of experimental relevance in the context of single molecule devices^ and of quantum dots^ and is important as a 
' paradigm problem for the development of techniques and insights. Two crucial issues in nonequilibrium physics are 
^T) ■ dephasing and decoherence. In the model we study dephasing arises because the wave functions in the two reservoirs 
CS| ' evolve in time at rates which differ by A/i whereas decoherence arises from the flow of energy and particles across 
ly-j . the system. An important issue in nonequilibrium physics is to develop methods which allow these effects to be 
' systematically analyzed. 

C — ■ ,, On the formal level, investigation of equilibrium systems is based on the partition function. Powerful techniques, 
most importantly the renormalization group method, enable one to eliminate putatively unimportant degrees of 
^ [ freedom and derive an effective theory governing the low energy behavior of interest. The renormalization group has 
been implemented in two intimately related ways: by considering changes in the self energies, vertex functions, and 
correlation functions in diagrammatic calculations, and by working directly with the partition function, generating 
an effective action describing only the degrees of freedom of interest. 

Most applications of renormalization group ideas to nonequilibrium problems have been based on the first approach: 
a diagrammatics is constructed using the Keldysh technique and then the flow of vertices, self energies and response 
functions under changes in cutoff is studied. In pioneering work Rosch and co-workers^ constructed a nonequilibrium 
scaling theory for the Kondo problem by identifying logarithms in perturbative calculations. Measurable quantities are 
computed, scaling equations for the coupling constants of the Hamiltonian were inferred from logarithmic dependences 
of observables on the upper frequency cutoff. In a very recent paper Borda et al^ used similar techniques to study 
the nonequilibrium behavior of the tunneling center model by perturbation theory in the dot-lead coupling. Paaske 
et. al provided further insigh t^°d^ into the physically crucial issue of decoherence rates, showing that voltage- induced 
decoherence enters differently into different observables, so that the analogy between temperature and decoherence is 
not precise. In a more recent set of approaches, a transformation is performed on the Hamiltonian itself. — 

While these approaches have established a number of basic results and concepts, the development of the subject 
remains incomplete. In equilibrium, defining a renormalization group transformation directly on the free energy 
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FIG. 1: Sketch of energy levels and cutoffs; band energy cutoffs chosen for each lead to be symmetrical about Fermi energy 
in that lead. (A) (upper panel): chemical potential difference A/i = fJ,L — fJ.R = 0. (B) (middle panel) Chemical potential 
difference Ajj, ^ but less than cutoff 2^~^. (C) (lower panel) Chemical potential difference A/^ greater than cutoff 2^^^; no 
dissipation processes possible. 



provided important insights into the meaning of the transformations and the formal structure of the theory. A similar 
analysis out of equilibrium should lead to valuable insights including a clearer understanding of decoherence and 
dissipation and a more precise definition of the charges in the renormalization group equations and the ability to 
map one problem onto another. In this paper we therefore examine a simple model, the tunneling center problem, 
from the effective action point of view. The "tunneling center" model is perhaps the simplest example of a wide class 
of quantum impurity models such as the Kondo problem or the spin-boson problem, which involve a finite number 
of local degrees of freedom coupled to reservoirs. The equilibrium behavior of these models is very rich, involving 
nontrivial correlated states, dynamically generated energy scales and quantum phase transitions The equilibrium 
physics revolves around a competition between formation of a quantum coherent state of the local degrees of freedom 
and the decoherence associated with coupling to the reservoir. We wish to characterize the effect of departures from 
equilibrium on this physics. 

Effective actions for nonequilibrium problems have been discussed by several authors; for a review see e.gJ^. Here we 
define and analyze a renormalization group transformation directly on the effective action. Our approach is somewhat 
similar to previous work of Konig and collaborators ; "'^^i who defined a renormalization group transformation directly 
on the equation of motion for the density matrix. These authors treated the dot-lead coupling via (self-consistently 
resummed) perturbation theory. We study a simpler model in which we are able to treat the dot-lead coupling exactly. 
Our formalism allows us to deal with decoherence and orthogonality physics on the same footing. We show that the 
dephasing scale A// defines a crossover beyond which both the physics and the formalism changes. We find that the 
effective theory valid for energy scales lower than A/i is characterized by a richer structure of charges than is the 
corresponding equilibrium theory. It also involves, as a new parameter, a decoherence rate. This arises physically from 
the flow of current through the system. It enters the theory as an additional parameter, which is subject to additive 
renormalization (as is the total energy in usual RG treatments) . This is discussed in more detail in section IIIIl 

The rest of the paper is organized as follows. In section |lT]we define the model we consider and obtain the effective 
action. In section|TTT]we derive the relevant scaling equation, and in section|lV]we present the solution and its physical 
content, among other things displaying the different physics of decoherence and dephasing. Section |V] is a summary 
and conclusion, which outlines implications for other problems. 



II. MODEL AND DERIVATION OF THE COULOMB GAS 



There are many realizations of the tunneling center modeli^. We consider a two-state system, which we represent 
in spin notation, linearly coupled to the density of electrons in an N-fold degenerate electronic level (creation operator 
djj with a — 1...N), which is itself hybridized with two (L and R) leads characterized by free-fermion statistics with 
possibly different chemical potentials i.e., A/i = fJ-L—f^'B^ / 0. We consider mainly temperature T = 0. We assume that 
the coupling to the leads preserves the "pseudospin" index a. The Hamiltonian may then be written as {Sz,x = \'^x,z 
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is a spin matrix) 

H = Hioc + Hiiath (1) 

Hioc = S,B + AtS^ + XDS, Y1 '^^o.dc (2) 



= 1..JV 



a=L,R,a=l..N'' a=L,B.,a=l..N 

(3) 

where 0<^^l = -| — 6'ii;< ^ and we have absorbed the mean hybridization strength into the variable e. The density 
of states per pseudospin a in the leads is D^^. For the two level system we choose the spin basis which diagonalizes 
the coupling to the d electrons, and parameterize this coupling by a dimensionless variable A and the lead density of 
states. The "magnetic field" B is the level splitting of the two level system and the parameter At gives the tunneling 
between the states. Note that if _B = the Hamiltonian has a particle-hole symmetry, so its energies are invariant 
under combined operations of changing the state of the two level system from "up" to "down" and changing particles 
to holes (c? <-> d^); this simplifies the algebra without changing the conclusions. 

The model requires an upper cutoff, for the energy integrals. It is most convenient to assume that in each 
lead the cutoff is symmetrical about the chemical potential in that lead, i.e. in lead a the energy integrals run from 
fia — to fia + (see Fig. [IJ. We assume that in the starting model the cutoff energies are much larger than the 
chemical potential difference or than the level splitting parameter B. 

Crucial parameters of the model are the nonequilibrium phase shifts introduced in^^ and defined by 



Sa — arctan 



A cos^ 9a 



1 - X ^ is.gn(/ia - ^a)Asin^ Oa 



(4) 



Note that the phase-shifts are not independent variables, but are related to each other as S^q — Sl + Sr = arctan A/ (1 — 
A^/4). Also note that while S^q is just the difference of the equilibrium phase shifts ±tan"^ ^ associated with the 
two states Sz = Sl,r are not the differences of the nonequilibrium phase shifts associated with the two states. 

The behavior of the two level system is specified by the reduced density matrix, given at time t in terms of an 
initial condition p(0) at time t = by 

PS = Tr^i [e-'^V(0)e^^*] (5) 

Here Tr^i indicates a trace over all of the electronic degrees of freedom. We follow Anderson, Yuval and Hamann^S 
and expand ps perturbatively in the "spin- flip" amplitude A^. The spin flip events are viewed as particles with 
"fugacity" /^[^At] and interaction determined by the trace over electrons. The new features are the need for two 
time contours and a dependence of the interaction on the chemical potential difference A^. 

A term in the expansion of the density matrix consists of n_ spin-flip events at times running from 0_ to t along 
the time-ordered contour, followed by n+ anti-time-ordered from t to 0+. Some examples are shown on the left side 
of Fig [21 the top panel has 7i+ = 4, n_ = whereas in the two lower panels rt+ = n_ = 2. Labeling the times in this 
"Keldysh order" (time ordered on the — contour and anti-time ordered on the -I- contour) we have for the diagonal 
components of the density matrix (cr = ±1 below represents the state of the impurity spin), 

{^\Ps{t)\a) = E„_+„,=e.e„(-0"-*"+(^k"-p5or"^k) (6) 



rt dtn_ f-tn_ dtn_-l r-t^ r-z u,[,„_+i u.Ln_-i^-2 |. l 

Jo Jo I ■■■Jo t\ |_jo I Jo I •■■■Jo 

(ATe)"-+"+exp[/o({t,})] 

The interaction /q is obtained by evaluating the Tr^i, in other words by solving the Keldysh problem of electrons in 
the time dependent potential specified by the spin fiips. In principle this is a multiparticle interaction depending on 
all of the times {ifc}. In the equilibrium problem an essentially complete analytical solution exists^S, showing that the 
interaction is pairwise, with a logarithmic time dependence and coefficients given by the product of the sign of the 
charges (i.e. whether the spin flips from up to down) and the changes in scattering phase shifts. In the nonequilibrium 
case general analytical expressions are not known. The available evidence, including solutions at times short and long 
compared to A/i, perturbative calculations^ and numerics^ suggests the following structure, which is slightly more 
involved than in equilibrium. 



pt dt„_+i ptre_+l tit„_-).2 ntn 
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FIG. 2: Examples of spin-flip events in the Keldysh two axis representation (left side), and in the single axis representation 
employing quantum and classical fields (right side). The charges rii are indicated on the Figure. In the top and bottom 
examples the fields q = +1 for all instantons. In the middle example q = +1 for the left hand pair of instantons and = —1 for 
the right hand pair. Although the first and the third examples have the same quantum field configuration, but the different 
positions of the times on the Keldysh contour leads to different phase factors. 



To specify the structure, it is convenient to collapse the two-contour problem onto a single time axis by defining 
classical {Vd = ^{Sz{t-) + Sz{t+)}) and quantum fields {Vq = ^{Sz{t_) - Sz{t+)}). We then have a four state 
system, in which either Vq — and Vd ~ ±XD/2 or Vd — and Vq = ±XD/2. Transitions between these states are 
instantons, which we label by an integer n = ±1 giving the sign of the change in the quantum field (this is just the 
usual Coulomb gas charge) and an integer q — ±1 corresponding to the sign of the quantum field in the region where 
it is non-zero. The right hand panels of Fig [5] show some examples. 

There are in principle 16 pairwise interactions between the 4 kinds of instantons, but calculations reveal a simpler 
structure (see Appendix A). We find that the interaction between instantons at times U and tj is pairwise, as in 
equilibrium. As in equilibrium the sign of the interaction between instantons i and j is determined in the usual way 
by the product of the charges nirij, while the magnitude is determined by the product of the scattering phase shifts. 
However in the general nonequilibrium case the phase shifts are complex and take value 6i if = 1 and 6* if = — 1 
(see Appendix A). Finally, with each interaction is a phase factor. The interaction has a logarithmic time dependence. 



given approximately by In 



iaij\ti — tj\+(^ 



with phase factor = ±1 determined by the Keldysh time-ordering of U and 



tj, so that in the upper panel of Fig [3 023 = 1 whereas in the middle panel 023 — —1. We find it convenient to 
combine the phase factors into an over- all phase <&. The result is 



n_ +n_|_ 



(7) 



J2'n„q,v, exp Y.j<', n,njCo{qt, Qj, A^O + 



The prime symbol on the sum above is to keep track of constraints such as, two charges of the same sign cannot 
appear more than twice in sequence, and ti^ = for an expansion involving the diagonal component of the density 
matrix. 

The longi^ and short^°-^^ time limits of the interaction function Cq are known. At short times t < 1/A/i 



5l + 5r 



ln(^ 



(8) 
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independent of quantum fields. At long times t > we have 

Coiq = = Afit > 1; A^^ « l){t) - {Coiq = g - = iV f 4 + 4 ) (I^M^I) + Tneglil 



A/it > 1; A/i^ < l)(i) = Co(g = g : 



N 



Sl + Sr 
n 

, \Sr\ 



In- 



Am?, 



2 , ln(|A/lt|) +r„eg|t| 

m. 1 



(9) 



(10) 



From Eqns. [5] and 1101 note that the dynamics is characterized by an exponential time decay, which we use to define 
the decoherence rate T„eq which plays a fundamental role in the subsequent analysis: 



\S" — 5" I 

Tneg = -fneqA^l = NAfl ^ ^ ^' = NAfl^ ( ^A^ COS^ Ol Sm^ Ol + 



(11) 



The physics expressed by the interaction Cq in the A/i^ ^ 1 limit is as follows: for times less than the dephasing scale 
tdephasing = (A/Lt)~^ onc has the equilibrium result: the two level system interacts with one coherent combinations of 
the two leads (c^ = cos 6*^0!^ +sin6'7?c^); the other combination decouples. The coupling leads to the usual power law 
interaction with exponents given by the coherent phase shift Seq — 6l + Sr. Note that S^q is independent of quantum 
fields. For times longer than the dephasing scale, one has a richer structure. The model is effectively a two channel 
model with separate couplings to left and right leads. The interaction between instantons acquires a dependence on 
the quantum field. At times longer than the decoherence scale the interaction is cut off altogether. 

These limiting forms suggest the following decomposition of the interaction: 



Co{q, q, t ~ t', A/i, i) - Qo(g, q')hQ (^) + '^'^''^^ A/i(t - t')^ + ^ineqKeq {^—^,A^l{t - t'))(12) 



where the "charges" or phase shifts Qo,m obey the property 



Qo,M{q^q') = [Qo,Af(-g, -g')]* 
Qo.Miq, -q) ^ Qo,Mi-q, q) => ^"^[(^o,Af (9, ~q)] = 



(13) 
(14) 



The precise form of the functions ho m neq depend on the cutoff scheme, but for t/^ 3> 1 and A/<^ <gc 1 (see Appendix 




hoit/0 = ln{\t\/0 
1 



In- 



AmC 



A/it > 1 



implying that for a model with Afi^ <^ 1, 



-{A^itf, Afit < 1 

TT 

|A/it|, A/ii>l 



Qo(+,+) = A(% + % 

Qm{+,+) = -^SlSr 

go(+,-) = A^(^ + ^) 
SlS*r + 6r61 



(15) 
(16) 

(17) 

(18) 
(19) 



(20) 
(21) 
(22) 
(23) 
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The coefficient 7„eg is independent of the quantum fields and is given by 



Ineq = N 



2-n 



(24) 



The first term proportional to Qo in Eq [12] represents the effect of processes in which an electron emerges from one 
lead and is scattered back into the same lead; it is independent of A/i. The term proportional to (5j\/ represents the 
effect of processes in which an electron is transferred from one lead to another; it depends on A//. Finally, as will be 
seen, the last term expresses the decoherence; it vanishes if the coupling is only to one lead. 

It is also interesting to consider the expression for Cq in a model in which A/x^ > 1. This situation arises after 
rescaling. In this case the regime A/xi < 1 is not defined, and for all times t > we find. 



with (for the model with A/i^ > 1) Qm — and 



, A^< > 1) = go 



In 



1^1 



(25) 



tan 



- 1 A cos 



(26) 



In this limit there is no decay because the theory has no real process which allows nonconservation of energy (see 
Fig.IIl case C). 

Comparison of Eqs [51 and 1251 reveals an important point. If we apply the usual renormalization process of reducing 
bandwidth we must pass from the model which gives rise to Eq[5]and contains decoherence to the model which gives 
rise to Eq [23 No decoherence processes exist in this latter model, but the renormalization maps one model onto 
another model with the same physical content. We therefore conclude that one consequence of renormalization must 
be the generation of a decay rate, which appears as an extra parameter, additional to what is directly computed from 
the small bandwidth model. 



III. DERIVATION OF SCALING EQUATIONS 

In the formulation given in Eq. [3 the nonequilibrium two level system is seen to be a function of the dimensionless 
parameters A^^, A//^ and 5l,r- In this section we construct a renormalization group analysis by following the 
usual procedure of reducing the energy cutoff, i.e. increasing the time cutoff from ^ to ^' = (1 + A)^, integrating 
out the degrees of freedom in the eliminated interval and determining the consequences for the remaining degrees 
of freedom. These effects were considered for the equilibrium problem by^'': reducing the energy cutoff (increasing 
the time cutoff) leads to changes arising from the ^ dependence of Cq. Rescaling leads to the simple "engineering 
dimension" changes (A^, A/i) (At, A/x)(^/^') k, (At, A//)(1 — A). Finally, some kink-antikink pairs fall within the 
time interval between ^ and ^(1 + A) and must therefore be removed from the renormalized theory. The procedure for 
the nonequilibrium problem is similar to the equilibrium problem, in that in the starting formulation the minimum 
separation in time between tunneling events is ^. In the rescaled theory, spin flip events separated by time intervals 
less than ^' cannot explicitly appear, but their presence will lead to a renormalization of the interaction between the 
processes which do explicitly appear in the theory. For small A-r^A the sequence of close tunneling events which 
appears with any probability are the "close pairs" shown in Fig. [3| A close pair may lie between two spin-flip events 
that occur on the same Keldysh axis (example pairs A and B in Fig [3] that lie between spin-flip events ii,i2), or 
may lie between spin-flip events on opposite Keldysh axes (examples C and D that lie between spin-flip events ^4)- 
Physically, a close pair corresponds to a dipole, and leads to a screening of the interaction between other spin flip 
events that may lie on the same Keldysh axis or on different Keldysh axes. Note that we always consider close pairs 
that lie on the same axis. A close pair with one member on each contour cannot be considered as a dipole; its removal 
would change the spin state at the final time t or initial time t = 0. 

Mathematically, the calculation is easiest to perform in the collapsed single axis representation in terms of classical 
and quantum fields. As an example consider Fig. [4] The interaction between the charge at ti and the charge at 
can get renormalized by integrating out eight kinds of close pairs, four of which are nearest neighbors to charge at ti, 
while the other four are nearest neighbors to the charge at tk- We show the four close pairs that are nearest neighbors 
to ti in Fig [H Close pair I can occur in two ways which have the same quantum field configuration, but different 
classical field configuration. 
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FIG. 3; Scheme for integrating out close pairs. 



The change in the density matrix, on integrating out the four close pairs is (A^ — iA^^) 



„2ri+2 



= P 



•In 



[l + A| (2 * / + // + ///)] 



(27) 



Summing over all possible positions of the relative separation e between the two charges of the close pair, and their 
center of mass position r one obtains. 



U ^ .?(1+A) 
ti-i ^ 



^exp 



/// 



exp 



Co(+, +, e) + X! (-C'o(+, qk, \tk~T ■ 

k 

~Co(+, +, e) + X! '^^"'^ (Co(+, \tk~T ■ 

k 

Co(-, -,e) + ^niTife (-Co(-,gfc, |tfc - r 



||) + Co(+,gfe|t,-T+||)) 



Co(+,gfe|tfe-T + -|) 



|) + Co(-,9fe|t,-T+||)) 



(28) 
(29) 
(30) 



In the above three expressions, the first term in the argument of the exponent represents the self interaction of the 
close pair, and the second term represents the interaction of the close pair with all other charges. The latter may be 
Taylor expanded in e. The term I on Taylor expansion takes the form 



Ju-i C 4 \ k 



dr 



(31) 



The integrals over r and e can easily be performed. The result is of 0(A), which implies that the change to the density 
matrix in Eq. [27]is 0(A|,A) ^ 1, and therefore can be re-exponentiated. Finally, integrating out the four close pairs 
of Fig m leads to the following renormalization of the density matrix, 

gAA|,e-<^o<«' ninfch2Co(+,9fc,|tfe-t.|)+2Co(+,gfc,|tfc-ti-i|)+Co( + ,gfcJtfc-t, + i|)-Co( + ,gfcJtfc-t,|)-Co(-,gfcJti+i-tfc|)+Co(-,gfcJifc-t,|)] 

Note that the term e"'-'''^^'' is from the self-interaction of the close pair. Since this term arises at short times i ~ ^, it 
has no explicit quantum field dependence, and therefore the quantum field label has been dropped. 

Repeating the above computation for all possible positions of close pairs, one finds that the initial factor of 



e ^ ^ 
lowing way 



2t^ 1 1 -2t,- 



cancel among each other, while the function CQ{qi,qk, \tk — ti\) is renormalized in the fol- 



Co(9j,g/c, \tk - U\) -> Co{qi,qk, \tk - U\) 
-AA2, [6Coiq,, qk, \tk -U\)~ Co(-(7„ qk, \tk ~ U\) - Co(g„ -<Zfe, \tk - U\)] e-^«(|=i) 



(33) 



Note that Co(t/C = 1) = if A^^ << 1 but r„egC if A^i^ > 1. 
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FIG. 4: Close pairs that are nearest neighbors to charge at time U 
The second effect arises from the exphcit dependence of the interaction Co on the cutoff f . We write, for infinitesimal 



A, 

with 



Co{t,^,A^l)^Coil,AA^lt) + SCo (34) 



Adding Eqns [33] and [53] leads to the fundamental scaling equation, 

= -Aye '""^ ' [6Co((7, q ,t)- Co{-q, q ,t) - Co{q, -q , t)] + -jj-^ (36) 

The second term in Eq. [3S] may be computed from the fundamental Eq. IA13I While a precise general expression is 
not known, the limits are established. For A/i^ <C 1, Eqns. [51 and [TUI show that -j^f^ is independent of time, with 
coefficient Qo + Qm independent of quantum fields. This is just the equilibrium scaling. For A/i^ 3> 1, Eq. 1251 shows 
that is again independent of time, with some coefficient Qq derived from scaling and not, in general, given by 
Eq. [26] The regime A/i^ ^ 1 requires a more careful treatment. Differentiation of the perturbative results derived 
in Appendix [B] yields two contributions: a function gM = which expresses the dephasing crossover by "turning 
off' the Qm contribution to Co as Afj,S, increases through unity, and an additional contribution proportional to t 
for A/it > f which expresses decoherence by cutting off the interaction between instantons. The decoherence term 



is characterized by a coefficient gneq = Ltt 



Fig [5] shows gneq calculated from Appendix [B] within 



1 d^w, 

'"eg t din? 

perturbation theory to 0{X^) and for hard and soft cutoffs. (Note that the hard cutoff model gives rise to oscillations 
which we have neglected). 

In summary, our scaling theory must keep track of the changes in the function Co as the cutoff is changed, but 
consideration of the short and long time limits shows that the model has effectively six coupling constants: "Coulomb 
gas charges" Qo(+i +)i Qo(+i ~) expressing the part of logarithmic interaction between tunneling events involving 
the same lead. Coulomb gas charges Qm{+, +), Qm{+, — ) expressing the part of the logarithmic interaction which 
depends on coherence between leads, a tunneling amplitude Axq (which acquire a quantum field dependence labeled 
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FIG. 5: Crossover function Qneq (in units of ^Qm^I-i) describing additive renormalization of decay constant r„eq, for hard 
(solid line) and soft (dashed line) cutoffs. 



by q), and the decoherence rate Tneq- These are renormahzed according to 



dA 



T± 



= At± 1 - - [Qo{±, ±) + Qm{±, ±)5M(AA*a] 



d^neq ( ; 



dln^ 

d^neq (~t~! 



dln^ 
dQo,M{+, 



dln£^ 
dQo,M{+, ' 



dln^ 
dQo,M{—, 



'4 r„e9(+, +)A|+ + - (r„eg(+, +)A|+ - r„eg(+, -)a: 

-4 r„e9(-, -)A|_ + i (r„eg(-, -)A^_ - Tnegi-, +)A 



=,-Co(l) 



-gneq{AfJ.^) 



~ [3r„eg(+, -)A|_ + 3r„,,(-, +)A|+ - r„eg(+, +)A|+ - r„e,(-, -)A^T-] e"^"^'^ 

5r„eg(A/iC) 



Qo,m(+, +)A|+ + - (Qo,Af(+, +)A|+ - Qo,m(+, -)A|_) 



-Co(l) 



= - [3Qo,Af(+, -)A|_ + 3Qo,Af(-, +)Ai.+ - Qo,m(+, +)A|.+ - Oo,m(-, -)A^-] e 



-Co(l) 



dln^ 



Qo.Aii-, -)A|_ + - (Qo.Af -)A|„ - Qo,Af(-, +)A|+) 



'Co(l) 



(37) 
(38) 
(39) 

(40) 
(41) 

(42) 
(43) 



Note that in these expressions Co(l) = ^neq^ for Afi^ ^ 1 and = otherwise. The meaning of Eq. [351 ISH HOI 
is that as the renormahzed chemical potential passes through the scale A/i^ = 1 an additive contribution to Tneq is 
generated. 



IV. SOLUTION OF THE SCALING EQUATIONS 

We now discuss the solution of the scaling equations. Notice from Eq. [^D] - [23] that in equilibrium (A/i = 0), the 
charges Qo.m are independent of the quantum fields. Moreover out of equilibrium, any explicit dependence on the 
quantum fields appears to O(X^) for asymmetric couplings, and to 0(A*) for symmetric couplings. We will present 
results that are valid to O(A^), and hence ignore the explicit quantum field dependence of the charges. In this limit, 
as we shall show, the dominant effect of voltage is due to the decoherence term. 

The physically relevant starting point for the solution of the scaling equations is one where the chemical potential 
difference is small compared to the cutoff scale. In this limit Co(l),g„e<j ~ 0, and 5a/(A^^) « 1; thus the effective 
Coulomb gas charge Q = Qo + Qm (physically, the level couples to the coherent combination of the leads). The 
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scaling equations thus become 

dAn 



dQ 



dln£, 



(44) 

4gA^ (45) 



dlnS^ 



4r„egA^ (46) 



Eqs 1441 and 1451 are the usual equilibrium scaling equations and are solved as usual; from the solution the behavior of 
Qm and Tneq is computed. The equations leave the combination / = Ay — 5(^5^ — ln{Q/2)) invariant. If / < 
then the model scales toward At = (i.e. is localized) while if / > then it is on the delocalized side of the phase 
boundary and At increases. We discuss the two cases separately. 

If the initial conditions are such that the model is on the localized side of the equilibrium phase diagram, then 
scaling proceeds until the cutoff crosses through the dephasing scale A/i^ — 1. Beyond this point, changes occur. 
First, the leads decohere: qm — * so the term proportional to Qm drops out of the scaling equations and the 
charge becomes Qo- Depending on the sign of Qm (i-e. the relative sign of and Sl) this may either make the 
system more localized or more delocalized. Second, and of greater significance, the decoherence rate Tneq acquires a 
positive additive term, arising from the function gneq in the scaling equation. Third, the effective coupling becomes 
(A^)^ = A|,e~'^°(^) « A|,e~^"=''^. Thus the theory at the scale f > (A/i)~^ is characterized by a fugacity InA^^, a 
decoherence rate Tneq, a charge Qo, and by the scaling equations 



dA^ 



T 



d\n^ 

dTjieq 

dln^ 
dQa 
dln£^ 



= A' 



1 _ ^ 
2 



(48) 

4r„egA2,e-'^«(i) (49) 
4goA|e~^«(i) (50) 



We see that scaling proceeds in the usual way until the nonequilibrium scale ^ = 1 /Tneq is reached; beyond this point 
the factor e^*^"^^^ cuts off the scaling and we are left with a perturbative theory. 

A particularly important special case occurs if J = 0. In this case the equilibrium fixed point is Ay = 0, Q = 2, and 
if the initial value of A/z is sufficiently small, this fixed point is approached very closely, so that at the scale ^ = 1/ A/x 

At(0 = -— (51) 

Q « 2 (52) 
Qm Q*M (53) 

Scaling through the crossover region then drives Qm Q*mi changes the basic charge to Qo = 2 — Q*m, generates a 
Tneq = f \Q*m\ Afj,, and does not change At significantly. Scaling then proceeds until ^ — > l/r„eq ~ 2/(77 \Q*m\ Afi). 
We therefore see that the dephasing crossover typically shifts the system away from the critical point, and that 
decoherence then cuts off the scaling. The decoherence cutoff occurs very rapidly, unless Q*m << 1, meaning that 
one of the leads is much more weakly coupled than the other one. In this case a significant nonequilibrium scaling 
regime can exist. 

If the initial condition is on the delocalized side of the equilibrium phase diagram, then again we distinguish two 
cases, according to whether or not the model flows to strong coupling before A/j,^ = 1 or not. In the latter case, the 
treatment outlined above applies. In the former case, the Kondo or coherence scale is larger than the dissipation rate 
and a different treatment, beyond the scope of this paper, is needed. 



V. CONCLUSIONS 



In this paper we have expressed the nonequilibrium tunneling center model in terms of a Coulomb gas defined on 
the Keldysh contour. The nonequilibrium problem has a richer structure than the corresponding equilibrium problem. 
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In particular the effective low energy theory is shown to be a Coulomb gas characterized by two parameters (charge 
and sign of quantum field). Crucial ingredients of the resulting theory are the dephasing arising because the wave 
functions in the two leads precess at rates which differ by the chemical potential difference, and a decoherence arising 
from the dissipative processes again allowed when the model is driven out of equilibrium. We showed explicitly that 
the decoherence effects cut off the power law interaction between instantons which is found at T = in equilibrium. 

Further, we generalized the standard equilibrium scaling theory of the model to the nonequilibrium case. We found 
that scaling through the dephasing crossover A^^ 1 generates an additive renormalization to the decoherence rate. 
From this we conclude that the decoherence rate is a fundamental parameter of the nonequilibrium theory, which 
must be explicitly considered in a renormalization process. We further showed explicitly how the decoherence cuts off 
the renormalization group flow. 

A few words on the generalization of this approach to other quantum impurity models, such as the nonequilibrium 
Kondo modeliSi. The key difference between the model studied here and the nonequilibrium Kondo model is the term 
AtSx in Eq. [2] responsible for spin-flip processes. The analog of AtSx for the Kondo model is {JLL + JRR){S+S-+h.c.) 
where S is the impurity spin, while s = '^^^^f^y V'lo-^V'fc'CT' are the electron spins which have been written as the 
following linear combination of the two leads ibka = , \ =CkaL H — 7^^^=CkcrR- Thus in the Anderson- Yuval- 
Hamann procedure applied to the Kondo model, the rt_ + n-i_-th order expansion in the spin flip amplitude involves 
the computation of 



Tr 



electron 



= e'='°(*^'*-*"-+"+)i(ti, i2 . (54) 
rather than the quantity needed for the model studied in this paper 



Tr, 



electron 



In equilibrium, the quantity L in Eq. 1541 (referred to in the literature as the open line contribution)^^ acquires the same 
structure as that of the closed loop part e'-^", namely that of a Cauchy determinant. Thus Egns. [5H and Eq . [551 and 
therefore the Kondo model and the tunneling center model may be related by a simple redefinition of the phase shifts. 
Out of equilibrium this analysis breaks down because the dephasing between the leads occuring at A^t > 1 gives 
rise to a two channel structure similar to that analyzed by Fabrizio et al3i and Vladar et alj^. A direct numerical 
evaluation or a mapping to an explicit two channel Kondo model (with decoherence) could be employed. 

This work suggests several generalizations. First, the long-time exponential cutoff suggests that a numerical es- 
timation of the perturbation series may be possible. Second, the key issue in seeing a wide nonequilibrium scaling 
range is to get the decoherence time to be very large compared to the dephasing time. This does not occur naturally 
in the simple two lead models we have studied. A search for models, involving for example three leads, where this 
separation of scales occurs more naturally, may be of interest. 

Wide classes of models have been studied in equilibrium by Hubbard-Stratonovich methods, in which the partition 
function is expressed as a sum over configurations of auxiliary fields. In the strong coupling limit of many quantum 
impurity models a small number of auxiliary field configurations are relevant and the physics is controlled by tunneling 
between themj2£ Generalizing this analysis to the nonequilibrium situation is an important open problem22., for which 
the methods developed here may be useful. A useful first step might be a comparison to the Bethe-ansatz solvable 
interacting resonant level model^l. 

Acknowledgments This work was supported bvl nsf-dmf 04313501 

APPENDIX A: DERIVATION OF EQ. \6\ 

We start from the Hamiltonian in Eq. [T] which we write as a sum of two parts: 

H = Ho + Hi (Al) 



de-cl^^c,ac + ]l - / de (cos 9acl^^d^ + h.c.) 

,„ ,„ * Q=L,_R,a=l..Ar"'Ma-l/4 



a=l..N a=L,R,a=l..N 



Hi = AtS, (A2) 

When At- = 0, the Hamiltonian is exactly solvable and represents a noninteracting resonant level hybridized with the 
leads. The Anderson- Yuval-Hamann approach involves a perturbative expansion in At, treating Hq exactly. This 
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procedure was originally carried out for the partition function; we apply it here to the time-dependent density matrix 
p{t) determined from an initial condition po(t) via 



p{t) = e-'^Voe'^* 
The reduced density matrix for the impurity spin is defined as 

Ps{t) = Treip{t) 



(A3) 
(A4) 



and is a 2 X 2 matrix whose diagonal elements give the probability of the spin S to be up or down, while the off 
diagonal elements contain information about phase coherence. Rewriting 



^-iHt ^ ^-iHotrp^-iJ*dtV(t) 

V{t) = e'""*Hie-'"'>* = At [e'"''* S^e~'"°*] = ArS^it) 
a perturbative expansion in V{t) of Eg . I A4I mav be carried out, yielding 

Psit) = J2 (-»)"-*"MAtC)"-+"+ 



(A5) 
(A6) 



(A7) 







' r' dt„_+i r 


Jo e . 




Jo C ^0 



-+1 dt 



n_+2 



) • ■ • Sx{tn_+l)] 

We assume the initial density matrix 



Po 



/ T T 

_ ( PosPlead. 



V PosPleads 



(A8) 



where t, | represent the direction of the local impurity spin, while Pi^^ids represents the steady state distribution of the 
electrons when the local spin is oriented along t / i • The effect of the spin flip term would be to modify the diagonal 
components of ps from its initial value, and also to introduce off diagonal terms. The perturbative expansion for the 
diagonal component of ps is (note Hq appearing below implies Hq corresponding to Sz —1 / i), 

(T \ps\ T) = 



J2 (-*)"-*"+ (AtO"~+"^ 

n-,n+=0,2A... 



(A9) 



dtr, 



^ 



dU 



dh 




' dtn_+i r 


Jo e . 




Jo C Jo 



dtn_+l /■*"-+^ dt„_+2 dtn^ + 



Tr. 



+ (-i)"-z"+ (AtO"""^"+ 



n_ ,n4, — 1,3,5.. 





- dtn_-l 


dh 




' dtn_+i /■*' 




dtn++n_' 






Jo i Jo 




"Jo 




.JO C Jo 


C '"Jo 


c J 


















+ 1 



The above may be written as 

(T \ps\ T) 

Jo Jo 



C ■ • ■ Jo 5 Jo c Jo 



(AlO) 



Tr, 



Pso' 

.(-^)"-^"+ (AtC) 



/o f /o 



n_ ,n_|_— 1,3,5. 
rt2 dti 1 r j-t dtn_+i r-t 
•■■Jo £ Jo £ Jo 



-1 dtn_ 



■fo 



Tr, 
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where 



Vy^it) = XD 



E diWdcit) (All) 

a=l...N 



Sz{t) in lAlll is the expectation value of the operator Sz and switches between {+/—)^ at times ti,i2 ■ • 

The first term in Eq. lAlOl {n± = 0, 2, 4 . . .) represents "out-scattering" , the second term {n± — 1, 3, 5 . . .) represents 
"in-scattering" . To evaluate the Tr^i in the out-scattering terms we write the lead states in the basis of scattering 
states appropriate to the static potential {XD/2)J2a=i n '^a{^)'^a{t). The potential V^{t) then alternates between 
the values —XD X]a=i n ^^'^ 0- Similarly to evaluate the Trei in the in-scattering term we write the lead 

states in the basis of scattering states appropriate to the static potential (—AD/2) J2a=i n '^a {t)da {t) so Vi alternates 
between XD J2a=i n da{t)da{t) and 0. The density of states of the two scattering problems is identical. 

Since, the leads are noninteracting electrons, with nonequilibrium imposed by Ha 7^ Mb, the evaluation of the Tr^i 
for a given configuration of spin-flips reduces to a problem of single particle quantum mechanics in a time dependent 
potential. We briefly outline the solution based on the nonequilibrium linked cluster theorem^^ which implies 



where 



E 



f dhTr[<l>g,g{h){G^s{h,h) + G^J{h,t,)}] (A13) 
Jo 9 Jo 



where for in-scattering, G^'^^{ti,ti) are Greens functions appropriate to the classical fleld (jid = 
^ iSz{t+) + S'^(t_) + 1), and the quantum fleld (j)g^g = ^ {Sz{t-) - Sz{t+)). They obey the Dyson equation 

=g + gi<j)cil + (bg,,T.,)G3 (A14) 



^G^ G^\ fq^ q^^ 

where G^ — ( „^ „\ ) ' -9 ~ ( n a]- Here g^'^^^ are the standard retarded, advanced and Keldysh Greens 



GfJ' ^ \0 g' 

functions of H, Eq[Tl with Sz = ^1/2- Note that g^'^ are short ranged in time and may be approximated as delta 
functions 

It is convenient to recast Ea lA14l as 

G' ^ 9 + gc/^g.gT^G' (A15) 



with g — (1 — 4>cig) ^ g- Explicitly, we have 



-R,Au (0d(O-AD/2)±z£) 



g^(t,t') = g"h-hg^ = -- ■ (A17) 

^ ^' ^ ^ {tD + XD/2~Mt))i-^D + XD/2-Mt')) ^ - 1' ^ ' 

with h ^ l/{t — t') the usual distribution function. Note that for times t,t' such that 0^ ^ 0, 4)ci — XD/2. 
Rearranging Eg lAlSl exphcitlv we find that G^ and obey the equations 



^99 



Eg lAlSl and Eq |A19l are singular integral equations. Noting that G^{t, t') ^ only when (t)q{t) ^ and that g^^^ are 
effectively delta functions, we see that the long time behaviors of G^ and G^ are the same. In equilibrium we may 
set /ii^jj — and write Eq I A19I explicitlv for i,t' separated widely in time. The important term is the last one, which 
is 



1 A 



d^i— X^G^(ii,t') (A20) 

t '-1 1 — 4- 
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where the prime denotes an integration only over those times ti for which 4>q{ti) ^ 0. From the standard properties 
of singular integral equations2fi'^ we identify the coefficient as the tangent of the phase shift, obtaining 



tan(5eg = 



A 



1 - 



(A21) 



and recovering the usual equilibrium Coulomb gas. 

In the non-equilibrium long time limit we follow Ng^^ and write — e"*^^(*^* )QK_s^^-nj.R(t-t )g^. We substitute 
this expression into Eq |A19[ write separate equations for Gf^ and G^, use Eq |A17l and note that for A/z(t — t') » 1 
the cross term between the e^^^^ term in and the e*'^^ term in gives an effective delta function contribution to 
the equation for G^ . This leads to a singular term of the form of Eq IA20I but with the phase shift replaced by 



i{m5a=L,F({t) 



(A22) 



Note from Eq. EH] that ,5a (-0gg) = -^aC'/'s?)- 

We briefly discuss the structure of the solution in the long-time nonequilibrium limit. From Eq. IA13[ we are 
eventually interested in the equal time limit of the Green's functions, which just as in equilibrium, have a divergent 
contribution due to the long time approximation made in deriving them. This issue can be resolved by solving for 
Gk.z by assuming that the Green's functions adiabatically follow the time dependent potential. The non-divergent 
part of Gk,z leads to the logarithmic interaction between the charges, which has the following form 



C/n 
- 
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a=L,R 



4>gq dSa 



TT d(l)g 



it) 



d\nXa{t) 
dt 



(A23) 



where 



Xa=L,R{t) = exp 



P 

2ni 



dt" , 1-itanSJt") 
In ^ ' 



t-t" l + tan^a(i") 

» Sajt") 



(A24) 



Using Eg. IA241 one may rewrite \nXa{t) ^ —^P J dt" jf^-^, further integrating Eq. IA23l bv parts one finds 



a=L,R 



it)^'-^ 

^ 't" ~t dt" 



(A25) 



Our model allows for a sequence of quantum fields which alternates between and ±^^^, and therefore may be 
written as, 



g\D 



q^[0{t-t2^-l)~e{t-t2^ 



(A26) 



i denoting times at which the quantum field changes, while, qi denoting the sign of the quantum field in the region 
where it is nonzero. From Eq. IA22I it follows that 



5a{t) 

dt 



J2 Saiq^)[Sit-t2^-l)~6{t~t2^)] 



(A27) 



In a similar manner 



Substituting Eq. [A27l and [A28l in Eq. lA25l we find 

^1 r QSa 













. d(j)gq 


i=l... 


. d(t>gq_ 



(qi) [0{t - t2^-l) - 0{t - t2i)] 



^'" = -i^/ t/^' ^ '^^^^ {q,mt~t2,-l)-B{t-t2,)]Y,5a{q^. 



^2i-l — t t2i 



(A28) 



(A29) 
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The integral over time t give rise to the logarithmic interaction between charges. Moreover the coefHcient of the 
logarithmic interaction between charges at ti and tj depend on the quantum fields Qi ^. In particular if Qi = qj = q, 
the coupling constant integral yields a coefficient 

On the other hand, if qi ~ —qj — q, the coefficient of the logarithm interaction between charges is 

a=L,R 

The overall ± signs before the coefficients of the logarithmic interaction essentially keep track of whether the spin has 
flipped up or down and may be used to define the Coulomb gas charge n.^. 

The discussion so far is valid for Afit ^ 1. In order to obtain an expression for Co for arbitrary Afit, we solve the 
Dyson equation perturbatively, and obtain expressions for Cq correct to second order in the scattering potential A. 
This is outlined in Section [B] The expression obtained interpolates the exact analytic expressions for A/it ^ 1 and 
Afit > 1. 



(A30) 



5 _ r r. 



iq) + Saiq) 



dS„ 



APPENDIX B: PERTURBATIVE EVALUATION OF Co FOR A SYMMETRIC AND HARD CUTOFF 



Let us turn to the evaluation of the time evolution operator at 0(A|n) (which corresponds to a single instanton in 
the quantum field), 



-Co(t) _ 



Pso 



(Bl) 



with Vj (f) defined in Eq. [XTTl 

By using the convenient representation da{t) = J^k a=L B.'^'^^^a.Ckaa{t), the expression for Co at O(A^) is 



Co{t)^^^ cos^OaCos^Op I dti I dt2 



al3=L,R 



^0 



*2(4„(il)c.,o(t2))(Ce./3(tl)c^,«(t2)) (B2) 



where < = f — ^i? < f ■ The limits of integration for ei_2 correspond to band edges that are abrupt and 
symmetrically located with respect to the chemical potentials (c.f. Fig[T]). At zero temperatures. 



Coit) 



NX 



n 



ti 



COS^ 9a COS^ dp / dti / dt2 

a0=L,R 



"'0 

t 



dei 



E 

al3=L.,R 



COS t/„ COS 



ti [i _ e-«(*i"*2)/«l 



^0 



(B3) 



It is also of interest to define a soft cutoff model with density of states pl,r{() ~ e ''^ mx-.h!?^ xhis gives Eg. IB3l but 



with 



fi-e-'<*i-*2)/ei 



Writing the above as a symmetric and anti-symmetric combination of the time arguments ti, t2 one obtains 

Co = i?e[Co]+/m[Co] 



(B4) 



where i?e[Co] is given by the symmetric combination of ti, <2 which after a straightforward evaluation of time integrals 
leads to, 



Re[Co]{t,Afj,,^) = NX ( 1 



2NX 



2NX- 



2 COS^ Ol COS^ 9fj 

cos^ 9l cos^ 9r 



[Fo(2t/0 - 2Fo{t/0] 



FoiAfit) 



FiiAfit) 



Fo(2t/^ + Afit) + F„{2t/^ - Afit) 



Fi(2t/f + Afit) + Fi{2t/C - Afit) 



- Fo{t/^ + Afit) - F„{t/^ - Afit) 



- Fiit/C + Afit) - Fi{t/^ - Afit) 
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with 



where Ci{z) = 7 + In z + dt- 



Fo{x) = -ln|a;| - 7 + Ci(a;) + (cos.T - 1) 
Fi{x) — xSi{x) 

Note Fq{x < 1) = -a;V2, Fi{x <Cl) = x^, Fi{x » 1) = ^x. 



The antisymmetric combination of ti,t2 leads to 

Im[Co{t)] oc it 



(B5) 
(B6) 



(B7) 



and represents unimportant energy renormahzation that vanishes for the particle-hole symmetric case. 

Identifying the coefficients above with perturbative expressions for the appropriate phase shifts defined in Eq. 31 
and defining the functions, 



Fo{2t/i + Afit) + Fo{2t/^ - Afit) 



FoiAfit) 



- Fo(i/e + A^it) - Fo{t/^ - Afit) 



(B8) 
(B9) 



Fi(A/xi) 



Fi(2t/f + Afj.t) + Fi{2t/^ - A^t) 



Fi{t/£, + Afit)~F^{t/^-Ant) 



(BIO) 



/iAf(i/e,A^i) 
one obtains Eq. [121 

To obtain the scaling function gneq giving rise to the additive renormahzation of r„eg, we differentiate Eq. lB3l in its 
soft cutoff analogue with respect to In^, and (to extract the long time behavior) t. The resulting integral may easily 
be evaluated. We plot the real part in Fig [51 
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